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ABSTRACT 

Light scattering by atmospheric dust particles is responsible for the polar- 
ization observed in some L dwarfs. Whether this polarization arises from an 
inhomogeneous distribution of dust across the disk or an oblate shape induced 
by rotation remains unclear. Here we argue that the latter case is plausible 
and, for many L dwarfs, the more likely one. Furthermore evolutionary mod- 
els of mature field L dwarfs predict surface gravities ranging from about 200 to 
2500 ms" 2 (corresponding to masses of ~ 15 to 70 Mj up i ter ). Yet comparison 
of observed spectra to available synthetic spectra often does not permit more 
precise determination of the surface gravity of individual field L dwarfs, lead- 
ing to important uncertainties in their properties. Since rot at ionally- induced 
non-sphericity, which gives rise to non-zero disk-integrated polarization, is more 
pronounced at lower gravities, polarization is a promising low gravity indicator. 
Here we combine a rigorous multiple scattering analysis with a self-consistent 
cloudy atmospheric model and observationally inferred rotational velocities and 
find that the observed optical polarization can be explained if the surface gravity 
of the polarized objects is about 300 ms~ 2 or less, potentially providing a new 
method for constraining L dwarf masses. 

Subject headings: brown dwarfs — stars: low mass — polarization — scattering 
— stars: atmospheres 
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INTRODUCTION 



Together with the lowest mass stars, brown dwarfs belong to the class of ultracool 
dwarfs whose relatively low atmospheric temperature and high pressure results in clouds 
of refractory compounds that in turn influence the entire atmosphere. The condensates 
(most abundantly iron and Mg-silicate grains) form near the atmospheric temperature ex- 
pected from chemical equilibri um as they efficiently depl e te the condensible species from 



the gas phase above the cloud (jAckerman fc Marley Il200ll ; lAllard et al. 1 1200 It iTsuji 112005 



Burrows et al. 112006c iHelling et al. II2008I ). As the objects cool over time, the dust eventually 
settles down gravitationally below the visible atmosphere. 

The observed spectra and photometry of ultracool dwarfs have been compared against 
predictions from models that incorporate the current unde rstanding of atmospheric physics, 



chem istry, dynamics and most importantly cloud processes (ICushing et al. Il2008t IStephens et al. 



2009). When combined with brown dwarf evolution tracks ( ISaumon fc Marley 



2008 



Baraffe et al. 



20021 ) these models can place fairly tight constraints on the effective temperature, T cS , of 
dwarfs with quality spectra. However, model fitting to date can leave the surface gravity 
of these objects poorly constrained, often by up to an order of magnitude. Since the radii 



of evolved brown dw arfs are only weakly dependent upon mass (IChabrier fc Baraffe 1 12000 



Burrows et al. 11200 it ), surface gravity is nearly directly proportional to mass for mature field 
L dwarfs older than several hundred million y ears. For such objects, sur face gravities lying in 
the range of 200 to 2500 ms -2 are expected ( ISaumon fc Marley II2008I ). As very few brown 
dwarfs have known dynamical masses, there has yet been few independent tests of the masses 
and temperatures derived from the spectral fitting. 

Given the state of both observational and theoretical constraints on dwarf surface grav- 
ity, an independent constraint on the surface gravity of ultracool dwarfs is sore ly needed- 
There are some spectral indicators of low mass, particularly for early type L dwarfs. ICruz et al. 
( 120091 ) identify weak alkaline absorption lines, and differing strengths of metal oxides and 
hydrides as compared to typical L dwarf spectra as signs of low surface gravity by analogy 
to spectra of giant stars. The physical processes underlying these unusual spectral features 
are as yet poorly understood and the value of surface gravity at which they become apparent 
is uncertain. Imaging polarimetry can provide another independent metric of constraint for 
surface gravity. 



Linear polarization, almos t certainly arising from dust scattering (jSengupta fc Krishan 



200 ll ; ISengupta fc Kwok M2005I ) . has been detected in the optical bands from a good number 



of L dwarfs covering almost the entire range of spectral types L0-L8 (IMenard et al. 



Zapatero Osorio et al. 1 120051 ; iTata et al. I |2009j). Observations by IMenard et al. I (120021 ) 



2002 



show that 25% of L0-L3 dwarfs and 50% of L3.5-L8 dwarfs in their sample are intrinsi- 
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cally polarized while iZapatero Osorio et al. I ( 120051 ) have detected polarization from 15 ±9% 
of L0-L3 dwarfs and 43 ± 17% of L3.5-L8 dwarfs in their samples. In principle the ob- 
served polarization could arise from the presence of magnetic field. However, radio, X-ray, 
ultra-vi olet and Ha observ ations point to a lack of magnetic activity in mature field L 
dwarfs (IBerger et al. 1120101 ); magnetic field strengths in the range 100-1000G have been 
deduced from observations of 8.3 GHz radio emission from a few brown dwarfs. This im- 
plies that synchrotron processes will not lead to significant linear polarization in the optical 
(IMenard et al. Il2002l ). Comparing the small net linear po larization (of order a few times 
0.01%) detecte d from a samp l e of A p stars (ILeroy I Il995l ) with about 1 kG dipolar field 
at the surface, IMenard et al. I (120021 ) pointed out that the observed optical polarization of 
ultracool dwarfs could not be explained by Zeeman splittin g of atoms or molecules. Fur- 
thermore the warmer M dwarfs are found to be unpolarized (IMenard fc Delfosse 1 120041 ) . M 
dwarfs have little or no atmospheric dust although they should have stronger magnetic field 
(IMohanty et al. 1120021 ). So, we conclude dust scattering polarization is the most plausible 
physical process that can account for the observed polarization of L dwarfs. 

Atmospheric dust can produce a net polarization if either the dust is spatially inhomoge- 
neous on large scales or the dust is homogeneous but the disk is oblate because of rotational 
distortion. We argue here that while the former mechanism cannot be ruled out, the latter is 
consistent with the current polarization observations if the polarized objects have relatively 
low gravity. A rigorous polarization survey of L dwarfs would distinguish between these 
mechanisms, providing either a new method to probe dust cloud morphology or 'weigh' the 
gravity of ultracool dwarfs. 



THE ATMOSPHERIC MODELS 



In order to test if plausible, spatially uniform dust clouds can reproduce the observed 
polarization of the L dwarfs, we employ a grid of one-dimensional atmosphere models 



(lAckerman fc Marlev Il200ll : Marlev et al. 112002 : iFreedman et al. 1120081 : ISaumon fc Marlev 



20081 ) for specified T e g and surface gravities g. The atmosphere model parametrizes the ef- 
ficiency of sedimentation of cloud particles through a scaling factor / se( j. For a fixed T e g, g 
and f se d the model uniquely predicts the variation in mean particle size and particle number 
density through the atmosphere which plays the crucial role in determining the scattering 
polarization. The atmosphere model employed here successfully reproduces the spectra and 
photometry of a large number of L dwarfs at a wide range of wavelengths covering near op- 
tical to mid-infrared regions as probed by ground and space-based telescopes. Model fitting 
typically constrains the effective temperature of an object of given spectral type within 100K 



-4- 



and generally rule out the case for f se <i = 1 (jCushing et al. 1120081 ; IStephens et al. 1120091 ) 



The gas and dust opacity, the temperature-pressure profile and the dust scattering 
asymmetry function averaged over each atmospheric pressure level derived by the atmo- 
spheric code are used in a multiple scattering polarization code that solves the radiative 
transfer equations in vector form to calculate the two Stokes parameter I and Q in a lo- 



cally plane-parallel me dium (ISengup 



Rayleigh phase matrix (ILiu fc Weng 



a fc Marlev 1 120091 ) . A combined Henyey-Greenstein- 



20061 ) is used to calculate the angular distribution of 



the photons before and after scattering. Finally, the angle dependent I and Q are integrated 
over the rotation-induced oblate disk of the object by using a spherical harmonic expansion 
method and the degree of polarization is taken as the ratio of the disk integrated polar- 
ized flux (Fq) to the disk integrat ed total flux (F r ). The det ail formalisms as well as the 
numerical methods are provided in ISengupta &: Marlev I ( 120091 ). 



ROTATION-INDUCED OBLATENESS 



We employ the Darwin- Radau relationship ( IBarnes fc Fortney 1120031 ) for estimating the 
rotation induced oblateness: 
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Here R e and R p are the equatorial and polar radii respectively, Q is the spin angular velocity 
of the object and K = I/(MR 2 )=0.261 for polytropic index n — 1 and 0.205 for n = 1.5, 
I being the moment of inertia. Comparisons with detailed structure models (D. Saumon, 
private communication) show that irrespective of their age, brown dwarf interiors can be 
adequately approximated by polytropes with 1 < n < 1.3 with the larger n being appropriate 
for higher gravities. As n increases, the oblateness decreases for a given rotational velocity 
and hence the degree of polarization decreases. In the present work we consider n = 1 and 
n — 1.5 as two extreme cases. 



4. RESULTS AND DISCUSSION 

4.1. Polarization from Oblate Dwarfs 

Using the foregoing modeling approach we computed disk integrated polarization for a 
variety of model assumptions. The effect of varying surface gravity, viewing or inclination 
or projection angle i, / se d, and rotation-induced oblateness on /-band (the bandpass at 
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which most of the data is available) polarization is presented in Fig.l for a fixed effective 
temperature T e g = 1800K. The degree of polarization p increases slowly with increasing 
oblateness and then increases rapidly for oblateness greater than about 0.18. This is because 
for relatively smaller oblateness the second harmonic in the spherical harmonic expansion 
is dominant, but as the oblateness increases the fourth and the sixth harmonics contribute 
significantly increasing polarized flux Fq. At the same time, higher harmonics become 
dominant and hence reduce the total flux Fj (more reddening due to limb darkening). As 
a result, p = Fq/Fj increases rapidly. However, as the inclination angle decreases, the 
variation of polarization with respect to the oblateness changes noticeably. All else being 
equal thicker clouds (smaller / sed ) produce greater polarization, especially when the surface 
gravity is high. 

Fig.l shows that the observed amount of L dwarf linear polarization can be produced by 
dust scattering only if the oblateness is greater than about 0.18 irrespective of any allowed 
value of the parameter^. In the absence of a dust cloud, polarization at /-band is negligible 
for any oblateness beca use Rayleigh scattering yie lds significant polarization only at shorter 



wavelengths (5-band) (ISengupta fc Marley 1 120091 ) . As shown in Fig. 2, the disk integrated 



degree of polarization p remains almost the same within the range of T e g 1800-1300K roughly 
corresponding to spectral types L3-L8 and falls rapidly at higher T e g, where clouds form at 
lower pressure and are thinner. At T e g < 1800 K, clouds are found deeper in the atmosphere, 
are optically thicker, and produce significant polarization. The transition from L to T dwarfs, 
i.e., from cloudy to cloudless atmosphere occurs at about 1300K as the clouds dissipate 
or settle below the photosphere. Above T e g ~ 2400 K there are few condensates. Thus 
polarization is a marker for the presence of substantial cloud layers. 



4.2. Comparison to Observations 



Menard et al. I (120021 ) detected confirmed polarization from three L dwarfs, marginal po- 



larization from two and no polarization from three L dwarfs. We consider the five confirmed 

J2005h 



Zapatero Osorio et al. 



and marginally polarized L dwarfs from this observation, 
found confirmed polarization in /-band from 9 L dwa rfs out of 33 target s. Out of these 9 L 
dwarfs, one (2MASS J1507-16) was also observed by iTata et al. I (120091 ) who detected con- 
firmed polarization from three L dwarfs. iTata et al. I (120091 ) detected polarizat ion as high 
as 5.2 ± 0.9 % in the /-band and 0.67 ± 0.17% in the #-band JTata et al. Il2009h of 2MASS 



1 The stability limit oblateness for uniformly rotating polytropes with n = 1.0 and 1.5 is 0.44 and 0.38 
respectively (James 1964). 
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J1731+27. This object also shows a very high H a equivalent width (-5.98) compared to the 
average value that is very small or zero. Spi tzer IRAC observations exclude a warm, but not 
a cold circumstellar disk (ITata et al. II2009I ). So, we exclude this object from consideration. 



For each remaining 15 objects, we computed the rotational velocity required to produce 
sufficient oblateness to reproduce the observed polarization for the T e g of the object (based 
on its spectral type). The spectral types, adopted T e g, photometric variability, detected 
/-band polarization along with the associated errors and the projected rotational velocity 
inferred from high resolution spectra as well as the same required to match the observed 
polarization are provided in Table 1. T e ff for almost a l l obje cts is derived from the optical 



spectral type by using equation (4) of Stephens et al. I ((20091). For DENIS-P J2252-17, the 



infrared spectral type is used as the optical spectral type is not known. T eS calculated from 
optical and infrared spectral type differs by less than 100K for all objects except for 2MASS 
J0141 + 18 which is a LI object in optical but L4.5 object in infrared and so its effective 
temperature r anges between 2100K a nd 1550K as derived by using equations (4) and (3) 



respectively of IStephens et al. I (120091 ). 



The projected rotational v elocity V sin(i) of a few L dwarfs showing confirmed polariza- 



tion i s inferred observationally (IMohanty &: Basri ll2003l : lReiners fc Basri 1120081 ; iBailer- Jones 



20041 ). However, the projection angle i is not known allowing a wide range of the values of V 
that along with g determines the oblateness. While a smaller value of i yields less polariza- 
tion, it gives rise to higher rotational velocity for a fixed ^sin(z) and hence more asymmetry. 
Taking the values of V sin(i) comparable to the observed values, we find that the models 
produce polarization comparable with that observed in /-band only when the surface gravity 
less than 1000 ms~ 2 and i < 45°, leading to substantial disk a symmetry. Fig.2A shows 
that all the five observational data points of iMenard et al. I (120021 ) can be fit well by setting 
i = 30° and g = 300 ms~ 2 with V sin (2) in range of 40 to 50 kms^ 1 which is within or slightly 
higher than the observed values. The observed /-band polarization of 2MASS J1507-16 fits 
well by using the observed V sin(i) = 27.2 km s -1 with g = 300 ms~ 2 and i = 30°. 

Except for Kelu-1 and 2 MASSW J1412+16, the r otatio nal velocity of most of the po- 
larized L dwarfs observed in IZapatero Osorio et al. I (120051) is unkn own. Fig. 2B shows 



that six out of the eight data from the IZapatero Osorio et al. I (120051 ) sample and one from 

2 MASS J1807+50 - whose obs erved projected rotational velocity is 



Tata et al. 


( 


2009) - 2 


76kms _1 ( 


Tata et al. 



20091 ; iReiners fc Basri 1 120081 1 . can be fit if V is as high as 90- 



105kms _1 and i = 90° at which the polarization is maximum. As this corresponds to 
a rotation period near an hour, a more likely explanation assuming homogeneous clouds 
would be an even lower surface gravity for these objects. For example, if = 80kms _1 , 
an equal amount of oblateness can be achieved by lowering the surface gravity to 240- 
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175 ms 2 . Alternatively, other physical processes such as surface banding or inhomogeneity 
m ay give rise to such high polari zation. We note however that the mean polarization value 



of IZapatero Osorio et al. I (120051 ) could be high because of the larger error bars owing to 



smaller telescope aperture than in the other studies (F. Menard, private communication). 

Among the polarized objects, Kelu-1 is most likely a triple system with uncertain po- 
larization contribution from the components. 2M ASS J2244+20 whi ch shows polarization of 



2.48 ± 0.47% is extremeley red both in the optical and near-infrared (IZapatero Osorio et al. 



20051 ). Hence the high polarization could be attributed to the presence of abnormally high 
amount of dust. The high polarization of 2MASS J1412+16, an L0.5 dwarf, while having 
small rotational velocity remains unexplained. Figure 2 shows that polarization does not al- 
ter drastically when the polytropic index is increased from n — 1 to n — 1.5. For n = 1.5, a 
slight increase in the rotational velocity is needed in order to fit the observed data. However, 
for the range 1.0 < n < 1.3, the observed data can be fit without altering the values of V, g 
or i. 

The rang e of T P ^ and the value of / se d adopted here overlaps with that derived from 
spectral fit by IStephens et al. I ( 120091 ) for the three common L dwarfs. All the data points in 
Fig. 2 can be fit with both / se d=2 and 3 because for g < 300 ms~ 2 , the polarization profile 
is not too sensitive to 3 > / se d > 2 as implied by Fig. 1. 

Monte Carlo simulation of the field substellar mass function indicates that objects in the 
12-75 Mj upiter mass range should greatly outnumber lighter objects in the solar neighborhood 
and objects below 12-13 Mj upi t er are expected to constitute a modest fraction of field L dwarfs 
(IBurgasser 1 120041 ) . Given the small sample sizes and lack of uniform selection criteria in the 
polarization surveys, it is premature to draw any conclusion from the relatively high fraction 
of low gravity objects w hich we find. Mean while the value of gravity below which the spectral 
indicators identified by I Cruz et al. I (120091 ) become apparent is not yet known, thus the lack 
of 'low gravity' spectral indicators in these objects is not necessarily indicative. 



4.3. Polarization from Surface Inhomogeneities 



Inhomogeneous distribution of atmospheric dust and Jupiter like bands (IMenard fc Delfosse 



20041 ) may also produce detectable polarization. Whether or not the clouds of L dwarfs 



are homogenous or patchy remains an open issue. Existing spectral models assume spa- 
tially uniform dust cl ouds and generally accurately reproduce observed L dwarf spectra 
(IStephens et al. 1120091 ) . Surface inhomogeneities can produce pho tometric variability and 40 



50% of L dwarfs are found to be variable without any periodicity (jGelino et al. 1120021 ; iKoen 
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2003 ). However, 2 MASS J22 24-01 2MA .SS J1108+68, 2MASS J1658+70 f lGelino et al. 
20021 ) and 2 MASS J1048+01 fkoen IboOa l are all variable L dwarfs with no detectable po- 
larization implying that inhomogeneities do not always p roduce significant po larization. On 



the other hand, 2M A SS J1412+16, 2MASS J0036+18 flGelino et al. 1 120021 ). and 2MASS 



J1507-16 (IKoen 1120031 ) are not variable implying a lack of large scale surface inhomogeneity 
(or a particularly favorable morphology and viewing geometry) but all of them are polarized. 
Another example is 2MASSW J1048+0 1 which is a variable ob ject with comparatively low 
projected velocity (Vsini = 17kms _1 ) (IReiners fc Basri II2008I ). If inhomogeneity produces 
detectable amounts of polarization then this object should show polarization irrespective of 
its low rotational velocity but it is unpolarized. 

An alternative hypothesis is that some dwarfs exhibit a uniform, banded appearance. 
Such objects might still be polarized but not be variable. We conclude that unlike the case for 
oblateness induced polarization-which naturally explains both the magnitude and variation 
with spectral type of polarization (Figure 2)-surface inhomogeneities require reliance on 
special cloud morphologies and viewing angles. Nevertheless both mechanisms likely play a 
role in producing polarization in some objects. 



5. CONCLUSIONS 

We have found that model L dwarf atmospheric structures which generally well repro- 
duce the spectra of known objects, predict full disk polarization of L dwarfs comparable 
to the values observed on some objects if the dwarfs are substantially oblate. Because the 
degree of oblateness varies inversely with gravity (Eq. 1), this mechanism requires fairly low 
surface gravities (g ~ 300 ms~ 2 ). Thus if oblateness is the primary mechanism by which L 
dwarfs become polarized, then polarization is a marker for low gravi ty. While some spe ctral 



indicators of low gravity have been identified in the literature (e.g. ICruz et al. 1 120091 ) . the 
actual gravity at which they become prominent has yet to be established. Thus polarization 
may serve as an indicator of moderately low L-dwarf gravity, at least for rapidly rotating 
dwarfs. 

A well constructed survey of a sample of field L dwarfs for polarization, V sin(z), variabil- 
ity, and spectral gravity indicators would test more definitively which mechanism (oblateness 
or surface inhomogeneities) is primarily responsible for L dwarf polarization. Once estab- 
lished, polarization would serve as a new constraint on the properties of newly discovered 
objects. However, this conclusion relies on the assumption that the adopted dust model 
describes the cloud distribution correctly. 
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Finally, we predict (Fig. 3) that the degree of polarization at J-band is comparable 
with that of /-band but is reduced at H- and K- bands. Again, in the infrared we expect a 
detectable amount of polarization only if the surface gravity is about 300 m s -2 or less. 
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T = 1800K 




0.1 0.15 0.2 0.25 0.3 

Oblateness 



Fig. 1. — Percentage degree of linear polarization calculated for the /-band as a function 
of oblateness. The numbers near the curves correspond to the value of the sedimentation 
efficiency parameter / se d- The polytropic index n = 1. 
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Fig. 2.— Model fits of the observed /-band polarization. The vertical error bars are 
observational errors and the horizontal ones are the spread of effective temperature for a 
particular spectral type. The numbers near the error bars correspond to the objects as 
listed in table 1. For all the cases / se d = 2. (A) The solid lines represent model with 
surface gravity g = 300 ms~ 2 and n — 1. From top to bottom they represent model with 
V = 48,41,48, and 27.2 km s -1 respectively. The dot-dash line represents model with 
n = 1.5 and V sin(i) = SOkms" 1 . The dash line represents that with g = 1000 ms -2 , n = 1 
and Vsin(i) = 48kms _1 . For all the cases except the one marked otherwise, % = 30°. (B) 
same as (A) but i = 90° and g = 300 ms -2 for all the cases. From top to bottom, the solid 
lines represent model with V sin(i) = 105, 100 and 90kms _1 respectively. The dash-dot line 
represents model with n — 1.5 and V sin(i) = 96 kms -1 . 
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Fig. 3. — Predicted J-, H-, and /f-band polarization of L dwarfs at different T eff and for 
different surface gravities. Only the model that fits the observed /-band polarization of a 
few L dwarfs is presented. Broken lines represent the polarization for /-band. 
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Table 1: Observed and derived quantities of polarized L dwarfs. 



Serial 


Discovery 


Spectral 


Estimated 


Photometric 


p% 


o % 


V sin(i) 


No. 


name 


type 


Teff(K) 


variability 






(kms -1 ) 






(Optical) 












1 


DENIS-P 


L8 


1340± 50 


Yes 


0.167 


± 0.04 


40.8± 8 a 




J0255-4700 




(1200-1300) 








(41-48) 


2 


LHS 102B 


L5 


1605 ± 100 


- 


0.105 


± 0.036 


32.5 a (41) 


3 


2MASS 


L3.5 


1780± 100 


No 


0.199 


± 0.028 


45 b 




J0036+1821 




(1700-1800) 








(48) 


4 


DENIS-P 


L2 


1966.5 ± 100 


- 


0.122 


± 0.042 


- 




J2036-1306 












(41-48) 


5 


DENIS-P 


L0 


2234 ± 100 


- 


0.083 


± 0.017 


- 




2000-7523 












(48) 


6 


2MASS 


L5 


1605.5 ± 100 


No 


0.036 


0.0 


27.2 c 




J1507-1627 




(1600-1700) 








(27.2) 


7 


2MASS 


LI. 5 


2100 ± 100 


- 


0.711 


± 0.142 


76 b 




J1807+5015 












(105) 


8 


Denis-P 


L7.5(IR) 


1419 ± 100 


- 


0.62 


± 0.16 


- 




J2252-1730 












(90) 


9 


2MASS 


L5 


1605 ± 100 


- 


0.58 


± 0.19 


- 




J0144-0716 












(90) 


10 


Kelu-1 


L2 


1966.5 ± 100 


Yes 


0.8 


± 0.27 


60 a 
















(105) 


11 


2MASS 


L1/L4.5 


1850 ± 250 




0.45 


± 0.15 






J0141+1804 












(90-105) 


12 


2MASS 


L0.5 


2166 ± 100 


No 


0.57 


± 0.19 


19 b 




J1412+1632 












B 


13 


2MASS 


L0.5 


2166 ± 100 




0.23 


± 0.06 






J1707+4301 












(105) 


14 


2MASS 


L4.0 


1718 ± 100 




1.38 


± 0.35 






J2158-1550 












(100) 


15 


2MASS 


L6.5 


1460 ± 100 


very red 


2.45 


± 0.47 





J2244+2043 (~ 105) 



Mohantv & Basri ( 


2003) 


'Reiners & Basri 


( 


2008) 


c Bailer- Jones ( 


2004) 



Note. 



Stephens et al. 



In colum n 4, the numbers inside brackets are T c g derived from synthetic spectra by 
(|2009h . p is the obse rved /-band mean po larization and a is the assoc iated error 



For objects 1-5, p and a are taken from iMenard et al. I ([20021 ). 6-7 from iTata et al. ] (l2009h . 8-15 from 



